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Abstract  This  work  examines  the  azimuthally  averaged  boundary  layer  structure  of  a  numerically  simu¬ 
lated  hurricane.  We  nominally  define  the  hurricane  boundary  layer  as  the  layer  in  which  the  effects  of  sur¬ 
face  friction  are  associated  with  significant  departures  from  gradient  wind  balance.  The  boundary  layer  in 
the  intensifying  primary  and  forming  secondary  eyewalls  is  found  to  be  nonlinear.  At  large  radii,  exterior  to 
the  eyewalls,  Ekman-like  balance  as  traditionally  defined,  is  found  to  hold  true.  Where  significant  departures 
from  Ekman-like  balance  are  found,  the  departures  are  characterized  by  large  vertical  advection  of  horizon¬ 
tal  velocity  through  the  depth  of  the  boundary  layer.  Shock-like  structures  are  not  found  to  be  prominent  in 
the  azimuthally  averaged  view  of  the  vortex  boundary  layer,  with  the  largest  azimuthally  averaged  radial 
gradients  of  the  radial  and  tangential  velocities  being  on  the  order  of  only  a  few  meters  per  second  per  kilo¬ 
meter.  Also,  in  the  radial  regions  of  the  eyewalls,  at  the  height  where  the  averaged  tangential  wind  is  a  max¬ 
imum,  the  radial  advection  of  radial  velocity  is  an  order  of  magnitude  smaller  than  the  agradient  force  per 
unit  mass.  Some  physical  implications  of  these  findings  are  discussed. 


1.  Introduction 

There  is  mounting  evidence  that  the  dynamics  of  the  hurricane  boundary  layer  is  a  key  for  understanding 
several  aspects  of  tropical  cyclone  intensification  and  mature  structure.  For  one  thing,  the  boundary  layer 
contributes  to  the  radial  distribution  of  moisture  and  convective  instability  via  sea-to-air  fluxes  of  latent  and 
sensible  heat  and  advective  transport.  For  another,  the  boundary  layer  dynamics  determines  not  only  the 
absolute  angular  momentum  of  the  air  leaving  the  layer,  but  determines  also  where  moist  boundary  layer 
air  erupts  out  of  the  layer. 

While  different  perspectives  could  be  adopted  to  define  the  hurricane  boundary  layer  [e.g.,  Zhang  et  a!., 
201 1],  at  an  elementary  level,  the  frictional  boundary  layer  of  a  tropical  cyclone  arises  on  account  of  the  fric¬ 
tional  stress  at  the  sea  surface.  In  this  work  "friction"  refers  to  the  turbulent  momentum  flux  divergence 
associated  with  the  near-surface  vertical  shear  and  the  surface  drag  stress  at  high  Reynolds  number.  In  com¬ 
parison  with  the  deep  vortex  flow  aloft,  the  boundary  layer  is  a  relatively  thin  transition  layer  that  connects 
the  bulk  vortex  with  the  frictional  surface. 

Of  course,  there  is  some  inherent  fuzziness  in  providing  a  precise  dynamical  definition  of  the  boundary  layer 
top.  The  ambiguity  arises  in  part  because  the  frictional  tendencies  are  likely  to  have  different  vertical  struc¬ 
tures  in  the  different  momentum  equations  and  also  because  the  agradient  force  (the  departure  from  gradi¬ 
ent  wind  balance)  is  not  strictly  zero  above  the  boundary  layer  so  defined.  Despite  such  fuzziness, 
dynamical  definitions  of  the  hurricane  boundary  layer  are  uncontroversial  in  the  outer  regions  of  a  hurri¬ 
cane,  where  convection  is  not  prevalent  and  there  is  subsidence  into  the  layer.  Flowever,  dynamical  defini¬ 
tions  of  the  boundary  layer  are  expected  to  have  limitations  in  the  inner-core  region  where  boundary  layer 
air  is  being  lofted  into  the  eyewall  clouds.  The  flow  in  this  region  is  akin  to  that  of  "separation"  in  aerody¬ 
namic  boundary  layers  [e.g.,  Anderson,  2005]  and  conventional  boundary  layer  theory  formally  breaks  down 
here  [Smith  and  Montgomery,  2010;  Montgomery  et  at.,  2014a].  It  is  difficult  to  unambiguously  define  the  top 
of  the  boundary  layer  in  this  region  because  significant  departures  from  gradient  wind  balance  are  still  pres¬ 
ent  that  are  not  directly  associated  with  the  surface  drag. 

In  this  study,  we  define  nominally  the  hurricane  boundary  layer  as  the  layer  in  which  the  effects  of  surface 
friction  are  associated  with  significant  departures  from  gradient  wind  balance.  More  specifically,  we  follow 
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common  practice  by  defining  the  boundary  layer  top  as  the  height  at  which  either  the  departure  from  gra¬ 
dient  wind  balance  vanishes  or  the  magnitude  of  friction  asymptotically  approaches  a  constant  (sometimes 
zero). 

Despite  the  widely  acknowledged  importance  of  the  boundary  layer,  there  are  differing  views  on  the 
dynamical  role  of  the  boundary  layer  in  the  spin-up  and  mature  structure  of  a  tropical  cyclone  vortex.  A  sur¬ 
vey  of  the  scientific  literature  reveals  an  eclectic  mix  of  ideas  being  debated  currently.  For  example,  in  an 
effort  to  explain  the  unprecedented  flight-level  observations  along  a  radial  flight  track  at  approximately 
500  m  altitude  in  the  rapidly  intensifying  Flurricane  Flugo  [Marks  et  al.,  2008],  some  research  has  proposed 
that  "shock-like"  structures — akin  to  the  classical  Burgers  equation  [e.g.,  Whitham,  1974] — are  essential 
elements  in  the  boundary  layer  dynamics  of  the  inner  core  of  an  intensifying  and  mature  tropical  cyclone 
[ Williams  et  al.,  2013;  Williams,  2015;  Slocum  et  al.,  2014].  Other  work  argues  that  a  quasi-linear  generaliza¬ 
tion  of  Ekman  theory  suffices  for  obtaining  a  zero-order  understanding  of  the  role  of  the  boundary  layer  in 
the  dynamics  of  secondary  eyewall  formation  [ Kepert ,  2013]. 

Still  other  work  argues  that  the  boundary  layer  dynamics  of  the  rapidly  rotating  region  of  a  tropical 
cyclone  are  intrinsically  nonlinear  [e.g.,  Smith,  1968;  McWilliams,  1971;  Smith  and  Montgomery,  2010; 
Montgomery  et  al.,  2014b]  and  cannot  be  captured  consistently  by  axisymmetric  balance  theory  [e.g.,  Bui 
et  al.,  2009;  Abarca  and  Montgomery,  2013,  2014]  or  quasi-linear  Ekman  theory.  In  this  nonlinear  view,  the 
boundary  layer  takes  on  a  new  dynamical  role  during  vortex  intensification  in  association  with  the  con¬ 
vergence  of  absolute  angular  momentum  ( M ).  Although  M  is  not  materially  conserved  in  the  boundary 
layer,  large  tangential  wind  speeds  can  be  achieved  there  if  the  radial  inflow  is  sufficiently  large  to  bring 
the  air  parcels  to  small  radii  with  minimal  loss  of  M.  This  spin-up  mechanism,  while  coupled  to  the  inte¬ 
rior  flow  dynamics  via  the  radial  pressure  gradient  at  the  top  of  the  boundary  layer,  is  tied  fundamentally 
to  the  dynamics  of  the  boundary  layer,  where  the  flow  is  not  in  gradient  wind  balance  over  a  substantial 
radial  span.  Indeed,  it  has  been  shown  using  theoretical  analyses  and  observational  data  that  this  bound¬ 
ary  layer  spin-up  mechanism  accounts  for  the  occurrence  of  the  maximum  tangential  wind  in  the  bound¬ 
ary  layer  [Bui  et  al.,  2009;  Sanger  et  al.,  2014;  Montgomery  et  al.,  2014a],  a  feature  that  has  been  found 
also  in  prior  observational  studies  [Montgomery  et  al.,  2006b;  Kepert,  2006;  Bell  and  Montgomery,  2008; 
Zhang  et  al.,  201 1]. 

The  differing  views  of  the  boundary  layer  summarized  above  motivate  basic  science  questions  about  the 
nature  of  the  dynamics  of  the  hurricane  boundary  layer.  For  example,  what  is  the  most  suitable  way  to  think 
about  the  dynamics  of  the  hurricane  boundary  layer?  Is  the  azimuthally  averaged  boundary  layer  structure 
of  a  hurricane  vortex  adequately  described  by  Ekman-like  dynamics,  or  shock-like  dynamics?  The  aim  of  this 
paper  is  to  contribute  toward  answering  these  questions  using  a  reasonably  high-resolution  "full-physics" 
mesoscale  simulation  of  a  mature  hurricane  vortex. 

An  outline  of  the  remainder  of  the  paper  is  as  follows.  Section  2  reviews  the  physics  of  Ekman-type  balance 
in  tropical  cyclone  vortices.  Section  3  introduces  our  methodology  to  investigate  the  dynamical  structure  of 
the  hurricane  boundary  layer.  Section  4  presents  and  discusses  the  results.  Section  5  offers  our  conclusions 
and  recommendations. 

2.  Ekman-Type  Models  for  Hurricane-Like  Vortices  Revisited 

As  described  by  Gill  [1982],  Holton  [2004],  and  McWilliams  [2011],  the  classical  Ekman  balance  comprises  a 
balance  between  three  horizontal  forces:  Coriolis,  pressure  gradient,  and  friction.  The  boundary  layer  flow  is 
characterized  by  a  large  Reynolds  number  and  is  therefore  turbulent.  By  friction,  then,  we  mean  the  diver¬ 
gence  of  the  Reynolds-averaged  eddy  momentum  flux  caused  by  the  surface  stress.  According  to  the 
boundary  layer  approximation,  this  frictional  force  per  unit  mass  may  be  approximated  by  the  vertical  diver¬ 
gence  of  the  mean  vertical  flux  of  horizontal  eddy  momentum.  Assuming  horizontal  homogeneity  and 
assuming  a  K-theory  turbulence  closure  to  parameterize  the  turbulence  effects,  the  horizontal  force  balance 
is  expressed  in  Cartesian  coordinates  as  follows: 
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1  dp  d 
p  dy  dz 


9v\ 

dz) 


where  the  indicated  variables  have  their  usual  meaning  and  K  is  a  vertically  variable  momentum 
diffusivity. 

For  curvilinear  flows,  such  as  a  rapidly  rotating  vortex  overlying  a  solid  frictional  surface  or  ocean,  (nonlin¬ 
ear)  Ekman-type  models  may  be  postulated  by  relaxing  the  horizontal  homogeneity  assumption  and  includ¬ 
ing  versions  of  the  apparent  forces  that  arise  from  the  circular  geometry  of  the  flow.  In  a  cylindrical  polar 
coordinate  system  for  a  stationary  axisymmetric  vortex,  an  Ekman-type  model  can  be  formulated  as  follows. 
Let  17  denotes  the  azimuthally  averaged  swirling  tangential  velocity  at  the  top  of  the  boundary  layer, 
assumed  to  be  in  gradient  wind  balance  with  the  azimuthally  averaged  pressure  field  P: 


r 


1  dP 
p  dr 


where  p  is  the  azimuthal  average  density  of  moist  air  and  r  is  radius  from  the  center  of  circulation.  Now  let 
v  denotes  the  departure  of  the  total  tangential  velocity  from  the  averaged  gradient  wind  velocity  and  let 
the  radial  velocity  in  the  boundary  layer  be  denoted  as  u.  Then,  for  small  departures  from  this  gradient  bal¬ 
ance  configuration,  the  first-order  boundary  layer  momentum  balance  in  the  radial  direction  is  between  the 
residual  cyclogeostrophic  force  and  radial  friction: 


-v 


d [ 
dz 


du\ 

dz) 


In  this  view,  departures  from  gradient  wind  balance  (left-hand  side)  are  balanced  by  the  vertical  gradient  of 
the  Reynolds-averaged  shear  stress  in  the  radial  direction  (right-hand  side).  To  connect  this  simplified 
momentum  balance  with  upcoming  language  used  in  section  3,  we  will  refer  to  the  negative  of  the  left- 
hand  side  of  this  equation  as  the  "Ekman  Agradient  Force,"  defined  as: 


In  the  tangential  direction,  the  first-order  momentum  balance  is  between  the  radial  flux  of  mean  absolute 
vorticity  and  the  vertical  gradient  of  the  Reynold-averaged  stress: 


u 


d_ 

dz 


The  foregoing  Ekman-type  model  has  been  used  frequently  in  the  hurricane  literature  because  it  appears  to 
capture  in  some  respects  certain  structural  features  of  the  hurricane  boundary  layer  flow.  These  features 
include  the  radial  distribution  of  the  vertical  velocity  at  the  top  of  the  boundary  layer  [e.g.,  Eliassen,  1971;  Eli- 
assen  and  Lystad,  1977]  and  the  existence  of  the  tangential  wind  maximum  within  the  boundary  layer 
[Kepert,  2001;  Kepert  and  Wang,  2001].  This  Ekman-type  model  is  still  being  used  as  a  theoretical  framework 
for  understanding  the  key  role  of  the  boundary  layer  in,  for  example,  secondary  eyewall  formation  and  evo¬ 
lution  [ Kepert ,  201 3;  Kepert  and  Nolan,  2014]. 


For  a  steady  state  boundary  layer  flow,  the  foregoing  system  neglects  the  nonlinear  advection  of  radial 
velocity  as  well  as  the  radial  and  vertical  advection  of  perturbation  tangential  velocity.  If  we  extend  our 
boundary  layer  to  include  the  viscous  sublayer  where  the  air  velocity  must  match  that  of  the  relatively 
small  (effectively  zero)  ocean  current  velocity,  we  note  here  that  because  v  must  be  of  the  size  of  V  to  sat¬ 
isfy  the  surface  boundary  condition,  and  u  is  comparable  to  v  because  of  the  expected  Ekman  layer  spiral, 
there  is  a  leading-order  inconsistency  by  including  the  linear  advective  fluxes  (u V,  vV)  while  neglecting 
the  nonlinear  ones  ( uu ,  uv,  vv).  This  reasoning  is  complimentary  to  the  scale  analysis  and  the  calculations 
presented  in  Vogl  and  Smith  [2009]  for  a  semislip  lower  boundary  condition  for  the  eddy  stress  at  ane¬ 
mometer  level  that  is  more  widely  used  in  geophysical  vortex  dynamics  [e.g.,  Eliassen,  1971;  Rotunno  and 
Emanuel,  1987],  Indeed,  we  will  see  below  that  some  of  these  nonlinear  horizontal  advective  terms  are 
often  significant  over  a  relatively  sizeable  depth  of  the  boundary  layer  momentum  balance,  as  are  the  ver¬ 
tical  ones  (uw,  vw). 
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3.  Analysis  Methodology 

The  methodology  of  this  work  consists  of  examining  the  kinematics  and  dynamics  (through  the  radial  and 
tangential  momentum  budgets)  in  a  mature  storm,  simulated  with  the  Regional  Atmospheric  Modeling  Sys¬ 
tem  (RAMS)  [Pielke  et  al.,  1992;  Cotton  et  at.,  2003].  The  analyzed  storm  was  originally  introduced  by  Terwey 
and  Montgomery  [2008]  and  was  revisited  by  Terwey  et  at.  [2013],  Abarca  and  Montgomery  [2013,  2014],  and 
Montgomery  et  at.  [2014b].  A  full  description  of  the  simulation  and  its  evolution  is  presented  by  Terwey  and 
Montgomery  [2008].  Here  we  highlight  that  the  simulation  uses  a  surface  flux  parameterization  based  on 
the  Louis  [1979]  scheme  (a  scheme  recommend  by  Kepert  [2012]);  the  radiation  scheme  introduced  by 
Harrington  [1997];  a  subgrid  scale  turbulence  scheme  based  on  Smagorinsky  [1963],  with  the  modifications 
from  Liliy  [1962]  and  Hill  [1974];  and  the  seven-species  microphysical  scheme  by  Walko  et  al.  [1995],  with 
the  specification  described  in  Montgomery  et  al.  [2006a],  While  we  focus  our  analysis  here  on  the  presented 
RAMS  integration,  all  the  conclusions  of  this  study  are  found  to  be  consistent  with  recent  diagnoses  by  the 
first  author  using  a  realistic  integration  of  Hurricane  Igor  (201 0).  The  Igor  simulation  was  performed  by  Davis 
et  al.  [2008]  using  the  Advanced  Research  core  of  the  Weather  Research  and  Forecasting  (AHW)  model.  The 
simulation  was  produced  by  The  Hurricane  Group  in  the  Mesoscale  and  Microscale  Meteorology  Division  of 
the  National  Center  for  Atmospheric  Research  (NCAR)  and  will  be  described  in  detail  in  an  upcoming 
publication. 

The  RAMS  simulation  is  an  idealized  integration  on  an  /-plane  (15°N),  with  sea  surface  temperature  held 
constant  (28°C).  It  is  executed  in  three  domains,  with  the  two  nested  domains  featuring  two-way  interac¬ 
tions.  The  nested  domains  are  located  at  the  center  of  the  corresponding  parent  domain.  The  domains  sizes 
are  168,  170,  and  251  grid  points  with  grid  spacing  of  24,  6,  and  2  km,  for  the  largest,  middle,  and  smallest 
domains,  respectively.  The  integration  is  performed  with  30  vertical  levels  (with  vertical  grid  spacing  varying 
from  300  m  near  the  surface  to  1.8  km,  at  about  26  km  altitude).  The  initial  condition  is  characterized  by  a 
quiescent  environment  ( Jordan's  [1958]  sounding);  a  weak  mesoscale  cyclonic  vortex  (10  m  s_1  maximum 
tangential  wind  speed  at  4  km  height  and  75  km  radius)  at  the  center  of  the  three  domains.  The  initial  vor¬ 
tex  is  assumed  to  be  in  thermal  wind  balance  and  possesses  a  positive  axisymmetric  moisture  anomaly 
near  the  center  (water  vapor  mixing  ratio  anomaly  gradually  increased  to  1 .3  g  kg-1  at  vortex  center). 

The  simulation  is  integrated  for  216  h.  Following  Terwey  and  Montgomery  [2008]  and  Abarca  and  Montgomery 
[2013],  hour  156  is  relabeled  as  hour  0.  Further  details  of  the  RAMS  experimental  setup  and  prior  analyses  can 
be  found  in  the  collection  of  papers  by  Terwey  and  Montgomery  [2008],  Montgomery  et  al.  [2006a],  Terwey 
et  al.  [201 3],  and  Abarca  and  Montgomery  [201 3]. 

To  optimally  analyze  the  nearly  axisymmetric  vortex,  we  remap  the  RAMS  variables  from  their  original  repre¬ 
sentation  in  a  Cartesian  coordinate  system  into  a  cylindrical  polar  coordinate  system.  The  center  of  the  coor¬ 
dinate  system  is  chosen  to  be  the  location  where  the  differences  between  the  azimuthally  averaged  radially 
profiles  of  sea  level  pressure  and  individual  profiles  taken  to  the  north,  south,  east,  and  west  were  mini¬ 
mized  following  the  methodology  of  Cram  et  al.  [2007].  We  note  that  the  results  exhibit  negligible  differen¬ 
ces  when  computed  with  the  center  obtained  as  the  centroid  of  the  potential  vorticity  field  in  the  lowest 
7.3  km  of  the  domain  (not  shown).  The  vertical  grid  spacing  of  the  cylindrical  grid  is  the  same  as  with  the 
Cartesian  grid.  The  radial  and  azimuthal  grid  spacing  are  2  km  and  1  degree  of  azimuth,  respectively. 

The  radial  and  tangential  momentum  tendency  equations  are  analyzed  in  their  respective  material  forms 
using  a  traditional  "azimuthal  mean  flow"  plus  "azimuthal  eddy"  partitioning  as  follows: 


du 

dt 


+  w  - 


du 

dz 


v'  du' 
r  dX 


+  w' 


du' 

dz 


V1  v'2  , 

—  +  —  +fv- 
r  r 


1  dp 
p  dr 


1  dp' 
p  dr 


+  Tr 


(D 


dv  _dv  _  dv  ,  dv '  ,  dv'  uv  u'v'  1  dp' 

—  +u  —  +w  —  +u'  —+w'— — —  —  - - — —-fu+F;. 

dt  dr  dz  dr  dz  r  r  rp  dX 


(2) 


Here  (u,  v,  w)  are  the  radial,  tangential,  and  vertical  wind  velocities;  r,  X,  and  zare  the  radial,  tangential,  and 
vertical  coordinates;  f  is  time,  p  is  air  density,  p  is  pressure,  and  Fu  and  Fv  are  the  radial  and  tangential  com¬ 
ponents  of  the  radial  and  tangential  momentum  subgrid  scale  eddy  flux  divergence,  computed  as  a  residual 
from  equations  (1)  and  (2).  As  is  usual  in  this  type  of  flow  partitioning,  overbars  represent  azimuthal  aver¬ 
ages  on  constant  height  surfaces  and  the  primes  denote  perturbation  (eddy)  values,  computed  as 
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Tangential  Velocity 


Figure  1.  Radius-time  plot  of  the  RAMS  azimuthally  averaged  tangential  velocity  (m  s_1) 
at  148  and  787  m  height.  The  dashed  black  lines  indicate  hour  4  and  hour  18,  where  the 
rest  of  the  study  focuses. 


departures  from  the  azimuthal 
averages.  (The  same  caveat 
]  expressed  by  Persing  et  at.  [2013] 
concerning  the  physical  mean¬ 
ing  of  the  "mean"  and  "eddy" 
terms  for  highly  localized  struc¬ 
tures  in  azimuth — e.g.,  with  vor¬ 
tical  hot  towers — applies  here  as 
well.)  The  terms  of  equations  (1) 
and  (2)  are  computed  from 
instantaneous  fields  (A  compari¬ 
son  of  the  results  using  1  and  2  h 
averages,  based  on  6  min  fre¬ 
quency  output — not  shown — , 
indicates  that  the  differences  are 
negligible  in  the  context  of  the 
conclusions  of  this  study.),  except 
for  the  time  tendencies,  which 
are  obtained  as  6  h  integrals  of 
centered  time  differences  (using 
6  min  output).  As  noted  above, 
the  frictional  tendency  terms  are 
computed  via  residual.  As  an 
example,  the  tangential  subgrid 
tendency  term  in  equation  (2)  is 
of  tangential  momentum,  mean  abso- 
terms  and  the  apparent  forces  and  azi- 


obtained  as  the  difference  between  the  sum  of  the  local  time  change 
lute  vorticity  flux,  the  mean  vertical  advection,  the  corresponding  eddy 
muthal  pressure  gradient. 

In  this  fully  nonlinear  context,  in  the  radial  momentum  equation,  departures  from  gradient  wind  balance 
associated  with  the  reduction  of  the  tangential  velocity  component  by  tangential  friction  are  accompanied 
by  an  effective  force  that  generally  does  not  balance  the  radial  friction  force.  The  specific  force  is  called 
here  the  agradient  force  and  is  defined  following  Smith  et  at.  [2009]  as: 


Ar.  1  dp  v2 
AP—  —  +  —  +fv 

p  or  r 


(3) 


In  the  subsequent  analysis,  we  will  confine  our  attention  to  the  lowest  3  km  of  the  inner  core  of  the  vortex, 
particularly  the  intensifying  primary  eyewall,  the  forming  secondary  eyewall,  and  a  region  radially  outside 
the  forming  secondary  eyewall. 


4.  Results 

4.1.  Boundary  Layer  Kinematics 

Figure  1  shows  the  radius-time  diagrams  of  the  azimuthal  mean  tangential  velocity  at  the  heights  of  1 48  and 
787  m.  It  is  similar  to  Figure  2  in  Terwey  et  at.  [201 3]  and  Figure  1  in  Abarca  and  Montgomery  [201 3,  2014].  The 
figure  summarizes  the  canonical  eyewall  replacement  cycle  of  the  simulation.  From  the  azimuthally  averaged 
perspective,  the  primary  eyewall  is  intensifying  up  to  hour  15  after  which  time  its  azimuthally  averaged  tan¬ 
gential  winds  weaken  monotonically.  Radially  outside  of  the  primary  eyewall,  there  is  a  progressive  amplifica¬ 
tion  of  the  azimuthally  averaged  tangential  winds  that  extend  beyond  the  domain  shown  (200  km).  From 
about  hour  24,  in  the  radial  region  of  the  tangential  velocity  amplification,  the  secondary  wind  maximum  of 
the  secondary  eyewall  emerges.  Finally,  the  figure  shows  the  radial  contraction  of  the  secondary  wind  maxi¬ 
mum  that  completes  eyewall  replacement  cycle.  In  this  study,  we  focus  on  hours  4  and  18,  when  the  primary 
eyewall  was  intensifying  and  when  the  secondary  eyewall  was  forming,  respectively. 

Figure  2  shows  the  radial  profiles  of  the  azimuthally  averaged  radial,  tangential,  and  vertical  velocity  compo¬ 
nents,  within  the  boundary  layer  of  the  storm's  inner  core  at  four  vertical  levels  between  148  and  1 138  m 
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Figure  2.  Azimuthally  averaged  tangential  (blue),  radial  (red),  and  vertical  (green)  velocity  components  (m  s-1),  as  a  function  of  radius,  at 
(left)  hour  4  and  (right)  hour  18  at  the  heights  of  (top  row)  1 138  m,  (second  row)  787  m,  (third  row)  457  m,  and  (bottom  row)  148  m.  For 
display  clarity,  the  average  vertical  velocity  component  is  multiplied  by  a  factor  of  5. 


height.  The  times  chosen  correspond  to  the  intensification  of  the  primary  eyewall  (hour  4)  and  during  sec¬ 
ondary  eyewall  formation  (hour  1 8).  These  times  have  been  found  to  be  representative  of  the  intensification 
of  the  primary  eyewall  and  the  formation  of  the  secondary  eyewall  within  the  first  20  h  depicted  in  Figure  1 . 

From  Figure  2,  we  see  that  the  largest  value  (and  radial  gradients)  of  the  radial  velocity  occurs  at  148  m  height 
(the  lowest  level  considered),  while  the  largest  tangential  velocity  occurs  at  787  m  height.  The  largest  vertical 
velocity  occurs  at  the  highest  level  shown  (1138  m).  Such  a  distribution  of  the  maxima  of  the  velocity  compo¬ 
nents  is  found  to  hold  true  for  the  intensifying  primary  eyewall  (hour  4,  left  plots)  and  the  forming  secondary 
eyewall  (hour  1 8,  right  plots).  The  vertical  location  of  the  maximum  values  of  the  horizontal  wind  components 
is  consistent  with  the  general  understanding  of  the  horizontal  momentum  balances:  the  strong  inflow  velocity 
is  primarily  attributable  to  the  effects  of  friction  and  it  reaches  its  maximum  magnitude  at  the  lowest  levels, 
where  vertical  momentum  diffusion  and  the  agradient  force  are  the  largest.  The  tangential  velocity  reaches  its 
maximum  near  the  top  of  the  frictional  boundary  layer,  where  the  radial  advection  of  tangential  momentum 
surpasses  its  depletion.  The  vertical  structure  of  the  horizontal  momentum  budgets  is  discussed  below. 

Throughout  the  boundary  layer,  the  largest  azimuthally  averaged  radial  gradients  in  radial  and  tangential 
velocity  occur  in  the  primary  eyewall,  both  during  the  intensification  of  the  primary  eyewall  (hour  4,  Figure 
2,  left  plots)  and  during  secondary  eyewall  formation  (hour  18,  Figure  2,  right  plots).  The  largest  azimuthally 
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Figure  3.  Vertical  profiles  of  (a)  the  three  velocity  components  and  (b-e)  the  radial  and  tangential  velocity  tendency  equations  at  hour  4.  The  quantities  presented  correspond  to  10  km 
radial  averages  radially  outward  of  the  radius  of  maximum  winds  (32-42  km  radius).  In  the  momentum  budgets,  the  dotted  lines  correspond  to  the  components  of  the  material  deriva¬ 
tive  and  are  plotted  as  positive  on  the  left-hand  side  of  equations  (1)  and  (2),  the  solid  lines  correspond  to  the  different  forces  in  the  budget  and  are  plotted  as  positive  on  the  right- 
hand  side  of  equations  (1)  and  (2).  See  text  for  details. 

averaged  radial  gradients  in  Figure  2  have  values  on  the  order  of  2  m  s_1  km-1  and  those  during  second¬ 
ary  eyewall  formation  are  at  least  an  order  of  magnitude  smaller.  Such  values  are  in  contrast  to  the  local 
radial  change  of  about  60  m  s_1  in  about  1  km  highlighted  by  Williams  et  al.  [2013],  as  they  were 
observed  in  Flurricane  Hugo  (1989).  Although  smaller  magnitudes  of  the  radial  gradients  are  expected  in 
an  azimuthally  averaged  view  than  those  obtained  from  in  situ  observations  along  a  radial  flight  segment 
through  the  eyewall,  the  smallness  of  the  azimuthally  averaged  radial  gradients  in  the  radial  region  of  sec¬ 
ondary  eyewall  formation  (and  even  during  the  radial  contraction  of  the  mature  secondary  eyewall  that 
are  about  1  m  s”1  km-1,  not  shown)  does  not  support  the  idea  that  shock-like  structures  are  relevant  in 
the  azimuthally  averaged  view  of  secondary  eyewall  formation  and  evolution  [cf.  Williams  et  al.,  2013; 
Williams,  2015;  Slocum  et  al.,  2014]. 


4.2.  Boundary  Layer  Dynamics 

Figures  3-5  show  the  vertical  structure  of  the  three  velocity  components  as  well  as  the  vertical  structure  of 
the  terms  in  equations  (1)-(3)  during  the  primary  eyewall  intensification  (Figure  3)  and  during  the  formation 
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Figure  4.  As  Figure  3,  but  for  hour  18  and  the  quantities  correspond  to  10  km  radial  averages  in  the  region  of  largest  tangential  wind  tendency  during  secondary  eyewall  formation 
(1 10-120  km  radius).  See  text  for  details. 

of  the  secondary  eyewall  (Figures  4  and  5).  The  terms  involving  perturbation  pressure  are  not  shown  (they 
were  found  negligible  compared  to  the  rest  of  the  terms  and  would  appear  as  vertical  lines  at  the  zero 
value  in  the  presented  diagrams).  Figures  3-5  show  10  km  radial  averages  and  their  results  are  consistent 
with  those  obtained  using  single  radius  and  5,  15,  and  20  km  radial  averages  (not  shown).  For  convenience 
of  presentation,  Figure  3  shows  radial  averages  between  32  and  42  km.  This  averaging  interval  is  10  km 
radially  outward  from  the  radius  of  maximum  winds  and  is  chosen  to  highlight  the  processes  leading  to  the 
emergence  of  the  maximum  tangential  wind  in  the  storm.  Figure  4  presents  radial  averages  taken  between 
1 10  and  120  km  radius,  focusing  on  the  region  with  the  largest  spin-up  tendency  associated  during  the  for¬ 
mation  of  the  secondary  eyewall  [see  Abarca  and  Montgomery,  2014,  Figure  5j].  Finally,  in  order  to  include  a 
representative  picture  of  the  storm's  dynamical  boundary  layer  structure  outside  the  strongly  convective 
region  of  the  vortex,  Figure  5  shows  radial  averages  taken  between  190  and  200  km  radius.  In  the  subse¬ 
quent  discussion,  we  refer  to  the  radial  region  of  Figure  3  as  the  intensifying  primary  eyewall  region,  to  the 
radial  region  of  Figure  4  as  the  secondary  eyewall  formation  region,  and  to  the  radial  region  of  Figure  5  as  the 
region  at  large  radii,  respectively. 

Figures  3a,  4a,  and  5a  highlight  similarities  and  differences  in  the  velocity  field  in  the  boundary  layer 
dynamical  structure  between  the  intensifying  primary  eyewall,  the  secondary  eyewall  formation  region,  and 
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Figure  5.  As  Figures  3  and  4,  but  for  hour  18  and  the  quantities  correspond  to  10  km  radial  averages  in  a  region  within  the  inner  core  but  radially  outside  of  the  eyewalls  (190-200  km 
radius).  See  text  for  details. 


the  region  at  large  radii.  As  proposed  in  section  1,  the  boundary  layer  top  can  be  nominally  defined  as  the 
height  where  the  disruption  of  gradient  wind  balance  associated  with  the  effects  of  surface  friction  van¬ 
ishes.  According  to  this  definition,  the  boundary  layer  top  is  seen  to  reside  at  a  height  of  about  2.0  km  in 
the  intensifying  primary  eyewall  (Figure  3b),  at  about  2.5  km  in  the  secondary  eyewall  formation  region 
(Figure  4b),  and  close  to  1.0  km  in  the  region  at  large  radii  (Figure  5b).  An  alternative  way  to  define  the 
boundary  layer  top  is  the  height  at  which  friction  asymptotically  approaches  a  constant  (sometimes  zero), 
inclusive  of  both  radial  and  tangential  momentum  balances.  This  second  definition  gives  a  boundary  layer 
height  of  about  2.5-3,  1.6,  and  1.0  km  height  for  the  intensifying  primary  eyewall,  the  forming  secondary 
eyewall,  and  large  radii,  respectively.  The  two  definitions  of  the  dynamical  boundary  layer  heights  are 
broadly  consistent  with  each  other. 

In  the  radial  regions  of  the  intensifying  primary  eyewall  and  forming  secondary  eyewall,  the  tangential 
velocity  attains  its  maximum  at  z  ~  787  m  and  ~1.1  km  height,  respectively.  Based  on  the  foregoing  diag¬ 
noses  of  the  boundary  layer  top,  it  is  clear  that  the  tangential  velocity  maximum  occurs  within  the  boundary 
layer  (Figures  3a  and  4a).  In  contrast,  in  the  radial  region  of  large  radii  the  height  of  the  tangential  velocity 
maximum  (~  1 .5  km.  Figure  5a)  is  marginally  higher  than  the  top  of  the  boundary  layer  (~1.0  km  according 
to  either  definition  considered).  As  the  radius  increases,  the  tangential  velocity  maximum  is  also  less 
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pronounced  relative  to  its  nearby  vertical  value  in  the  bulk  vortex.  From  Figures  3a,  4a,  and  5a,  the  differ¬ 
ence  in  magnitude  of  the  tangential  velocity  between  the  local  maximum  and  the  value  found  at  3  km 
height  is  8,  6,  and  2  m  s-1,  respectively. 

In  the  three  radial  regions  highlighted  here,  the  radial  velocity  is  characterized  by  a  layer  of  relatively  strong 
inflow  in  the  lowest  levels.  This  inflow  reaches  its  largest  value  at  the  lowest  level  shown  (148  m).  The 
regions  of  the  intensifying  primary  eyewall  and  secondary  eyewall  formation  contain  a  layer  of  outflow 
immediately  surmounting  the  layer  of  strong  inflow.  Such  an  outflow  layer  is  more  pronounced  in  the  inten¬ 
sifying  primary  eyewall  region  (outflow  velocity  of  6  m  s-1)  than  in  the  region  of  the  forming  secondary 
eyewall  (outflow  velocity  of  2  m  s_1).  A  corresponding  outflow  layer  at  large  radii  is  not  present  at  all.  In  the 
three  regions,  the  vertical  velocity  increases  with  height,  up  to  a  value  of  ~2ms“'  in  Figure  3a,  ~0.4  m  s^1 
in  Figure  4a,  and  ~0.3  m  s-1  in  Figure  5a. 

Ekman-like  balance  is  found  at  large  radii  (Figures  5b-5e).  In  the  radial  momentum  equation  (Figures  5b 
and  5c),  there  is  an  approximate  balance  between  the  subgrid  scale  eddy  flux  divergence  of  radial  momen¬ 
tum  and  the  inward  directed  agradient  force,  i.e.,  the  Ekman  Agradient  Force  as  defined  in  section  2.  In  the 
tangential  momentum  equation  (Figures  5d  and  5e),  there  is  an  approximate  balance  between  the  subgrid 
scale  eddy  flux  divergence  of  tangential  momentum  and  the  radial  flux  of  mean  absolute  vertical  vorticity 
(which  is  the  result  of  adding  the  centrifugal  and  the  Coriolis  forces  and  the  radial  advection  of  tangential 
momentum  multiplied  by  -1).  In  this  radial  region,  the  mean  and  eddy  components  of  the  radial  and  verti¬ 
cal  advection  of  radial  momentum,  as  well  as  the  mean  and  eddy  components  of  the  vertical  advection  of 
tangential  momentum,  are  relatively  small.  Thus,  an  Ekman-like  balance,  like  the  one  originally  proposed  by 
Eliassen  [1971],  is  found  here. 

Unlike  the  large  radii  radial  region,  Ekman-like  balance  is  not  found  in  the  regions  of  the  intensifying  pri¬ 
mary  eyewall  or  the  forming  secondary  eyewall.  In  both  of  these  regions,  the  subgrid  scale  eddy  flux  diver¬ 
gence  of  radial  momentum  in  the  radial  momentum  equation  (Figures  3b  and  4b)  is  too  small  to  balance 
the  agradient  force  throughout  the  depth  of  the  boundary  layer.  The  departures  from  Ekman-like  balance 
are  particularly  large  at  the  heights  where  the  agradient  force  is  directed  radially  outward.  The  tangential 
velocity  is  supergradient  here  [ Sanger  et  at.,  2014;  Smith  et  al.,  2009]  and  the  maximum  tangential  velocity  is 
found  in  this  region. 

Similarly,  in  the  tangential  momentum  equation  (Figures  3d  and  4d),  the  subgrid  scale  eddy  flux  divergence 
of  tangential  momentum  is  too  small  to  balance  the  radial  flux  of  mean  absolute  vorticity  through  the 
depth  of  the  boundary  layer.  Such  significant  departures  from  Ekman-like  balance  found  in  the  intensifying 
primary  eyewall  and  the  forming  secondary  eyewall  are  characterized  by  large  mean  vertical  advection  of 
mean  horizontal  velocity  through  the  depth  of  the  boundary  layer  (Figures  3b,  3d,  4b,  and  4d). 

Our  finding  of  the  dynamical  importance  of  the  aforementioned  nonlinear  acceleration  terms  is  consistent 
with  the  scale  analysis,  calculations  and  physical  interpretation  provided  by  Vogl  and  Smith  [2009].  These 
findings  do  not  support  the  hypothesis  that  the  dynamics  of  the  boundary  layer  in  the  eyewalls  of  a  mature 
storm  is  adequately  captured  by  Ekman-like  theory  [Kepert,  2013;  cf.  Montgomery  et  al.,  2014a]. 

Another  point  to  examine  here  concerns  the  structure  of  the  mean  radial  advection  comprising  part  of  the 
radial  acceleration  near  the  height  of  tangential  wind  maximum.  In  the  radial  regions  of  the  eyewalls,  at  the 
height  where  the  tangential  wind  is  a  maximum  (see  Figures  3a  and  4a),  the  mean  radial  advection  of  mean 
radial  velocity  is  an  order  of  magnitude  smaller  than  the  agradient  force  (Figures  3b  and  4b).  This  finding 
does  not  support  the  idea  that  the  mean  radial  advection  of  mean  radial  velocity,  and  corresponding 
"shock-like"  structures,  are  essentially  related  to  the  generation  of  supergradient  winds  in  the  boundary 
layer. 

There  is  one  additional  and  noteworthy  finding  that  emerges  in  our  analysis  of  the  azimuthal  mean  bound¬ 
ary  layer  structure  of  a  mesoscale  simulation  of  a  mature  hurricane  vortex.  For  the  three  radial  regions  con¬ 
sidered,  the  radial  subgrid  scale  force  exceeds  that  of  the  corresponding  tangential  subgrid  scale  force  near 
the  lower  surface.  In  retrospect,  this  is  perhaps  not  surprising  given  the  large  vertical  shear  of  the  inflow 
velocity  near  the  lower  surface  in  comparison  to  the  corresponding  vertical  shear  of  the  tangential  velocity. 
However,  it  would  certainly  be  incorrect  to  omit  the  radial  subgrid  scale  force  in  favor  of  the  tangential 
subgrid  scale  force  in  theoretical  attempts  to  analyze  the  boundary  layer  dynamics  of  a  hurricane  vortex 
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[cf.  Stern  et  at.,  2015].  The  full  implication  of  the  finding  that  the  radial  frictional  force  is  comparable  to  or 
larger  than  the  tangential  frictional  force  near  the  lower  surface  awaits  further  analysis. 

In  closing  the  results  section  of  this  paper,  we  reiterate  that  while  we  focus  our  analysis  here  on  the  pre¬ 
sented  RAMS  integration,  all  the  conclusions  of  this  study  are  found  to  hold  true  in  an  analogous  analysis 
based  on  a  realistic  WRF  integration  of  Hurricane  Igor  (2010)  having  a  lowest  model  level  of  50  m.  This  com¬ 
panion  analysis  will  be  presented  in  due  course. 

5.  Conclusions 

In  this  paper,  we  have  examined  the  azimuthally  averaged  boundary  layer  structure  of  a  numerically  simu¬ 
lated  hurricane  that  underwent  secondary  eyewall  formation.  This  analysis  aims  toward  a  clarification  on 
the  dynamical  nature  of  the  hurricane  boundary  layer.  Examples  of  contrasting  views  on  the  dynamical  role 
of  the  boundary  layer  are:  "Shock-like"  structures — akin  to  the  classical  Burgers  equation;  quasi-linear  gener¬ 
alizations  of  Ekman  theory;  and  a  fully  nonlinear  conception.  This  study  focuses  on  three  key  radial  regions 
of  the  boundary  layer  of  the  simulated  mature  hurricane:  the  intensifying  primary  eyewall,  the  forming  sec¬ 
ondary  eyewall,  and  large  radii  (considered  at  about  200  km  radius). 

We  defined  nominally  the  hurricane  boundary  layer  as  the  layer  in  which  the  effects  of  surface  friction  are 
associated  with  significant  departures  from  axisymmetric  gradient  wind  balance.  We  specifically  considered 
two  plausible  definitions  for  the  boundary  layer  top.  The  first  definition  is  the  height  at  which  the  departure 
from  gradient  wind  balance  vanishes.  The  second  definition  is  the  height  at  which  the  magnitude  of  hori¬ 
zontal  friction  asymptotically  approaches  a  constant  (sometimes  zero).  These  dynamical  definitions  were 
found  to  give  comparable  boundary  layer  heights.  The  azimuthally  averaged  tangential  wind  maximum 
was  found  to  be  located  well  within  the  boundary  layer  of  both  primary  and  secondary  eyewalls  and  found 
to  be  located  marginally  above  the  boundary  layer  at  large  radii.  In  all  radial  regions  considered,  the  maxi¬ 
mum  radial  velocity  is  located  at  the  lowest  level  considered.  These  modeling  results  are  consistent  with 
the  observational  findings  of  Kepert  [2006],  Montgomery  et  at.  [2006b],  Sanger  et  at.  [2014],  and  Montgomery 
et  at.  [2014b]. 

Although  Ekman-like  balance  is  found  at  large  radii,  it  is  not  found  in  the  regions  of  the  intensifying  primary 
eyewall  and  the  forming  secondary  eyewall.  At  large  radii,  there  is  an  approximate  balance  between  the 
subgrid  scale  eddy  flux  divergence  of  horizontal  momentum  and  the  inward  directed  agradient  force  in  the 
radial  momentum  equation,  and  the  radial  flux  of  mean  absolute  vorticity  in  the  tangential  momentum 
equation.  In  the  regions  of  the  intensifying  primary  eyewall  and  the  forming  secondary  eyewall,  Ekman-like 
balance  is  not  found.  In  these  radial  regions,  the  subgrid  scale  eddy  flux  divergence  of  horizontal  momen¬ 
tum  is  too  small  to  balance,  respectively,  the  agradient  force  in  the  radial  momentum  equation,  and  the 
radial  flux  of  mean  absolute  vorticity  in  the  tangential  momentum  equation.  Significant  departures  from 
Ekman-like  balance  are  characterized  primarily  by  large  mean  vertical  advection  of  mean  horizontal  velocity 
through  the  depth  of  the  boundary  layer.  The  dynamical  importance  of  these  nonlinear  terms  is  consistent 
with  the  scale  analysis,  calculations,  and  physical  interpretation  provided  by  Vogt  and  Smith  [2009]  who  con¬ 
ducted  a  test  of  the  validity  of  the  Ekman-like  model  using  analytical  solutions  for  the  quasi-linear  Ekman 
layer.  In  the  full  physics  simulation  examined  here,  eddy  advection  terms  contribute  also  to  the  departures 
from  Ekman-like  balance,  albeit  to  a  lesser  degree  than  the  mean  advection  terms  (that  have  larger 
magnitude). 

Shock-like  structures  are  not  found  to  be  prominent  in  the  azimuthally  averaged  view  of  the  hurricane's 
boundary  layer.  The  largest  azimuthally  averaged  radial  gradients  in  radial  and  tangential  velocity  occur  in 
the  primary  eyewall,  and  are  about  2  m  s-1  km-1.  Also,  in  the  radial  regions  of  the  eyewalls,  at  the  height 
where  the  tangential  wind  is  a  maximum,  the  mean  radial  advection  of  mean  radial  velocity  is  an  order  of 
magnitude  smaller  than  the  agradient  force  per  unit  mass. 

Additionally,  we  documented  that  near  the  lower  surface,  the  radial  subgrid  scale  force  exceeds  that  of  the 
corresponding  tangential  subgrid  scale  force.  The  findings  of  the  present  study  highlight  the  intrinsically 
nonlinear  nature  of  the  boundary  layer  in  the  simulated  hurricane.  The  findings  do  not  support  the  pro¬ 
posed  ideas  that  the  fluid  dynamics  of  the  boundary  layer  in  the  eyewalls  of  the  storm  are  adequately  cap¬ 
tured  by  Ekman-like  theory  or  that  shock-like  structures  are  relevant  in  the  azimuthally  averaged  view  of 
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secondary  eyewall  formation.  Finally,  we  noted  that  while  we  focus  our  analysis  on  the  presented  RAMS 
integration,  all  the  conclusions  of  this  study  are  consistent  with  those  from  an  analogous  analysis  based  on 
a  realistic  WRF  integration  of  hurricane  Igor  (2010).  Such  a  companion  analysis  will  be  presented  in  due 
course. 
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